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Abstract. The aim of this study was to evaluate the 
association between thymidylate synthase (TS), methylenetet- 
rahydrofolate reductase (MTHFR) and reduced folate carrier 
(SLC19A1) gene polymorphisms and the treatment efficacy of 
pemetrexed-based chemotherapy in advanced non-small cell 
lung cancer (NSCLC). Advanced NSCLC patients received 
pemetrexed and cisplatin every three weeks. Polymorphisms 
in the TS, MTHFR and SLC19A1 genes were detected in 
peripheral blood samples using DNA sequencing and Taqman 
PCR. An analysis of gene polymorphisms was performed 
with respect to the progression-free survival (PFS), response 
rate (RR) and overall survival (OS) of patients treated with 
pemetrexed. The median PFS times for patients with the TS 
2R/2R, 2R/3C or 3C/3C genotypes were significantly longer 
than those of patients with the 2R/3G, 3C/3G or 3G/3G geno- 
types (P=0.036). Patients with the SLC19A1 CC genotype had 
a significantly longer median OS compared with individuals 
with the homozygous and heterozygous genotypes (12.2 vs. 8.9 
and 7.3 months, respectively; P=0.022). The PFS and OS did 
not differ for the three genotypes of MTHFR assessed. The 
RR was higher in patients with the TS 2R/2R, 2R/3C or 3C/3C 
genotypes than in the other groups (P=0.044). The polymor- 
phisms of the 5'-UTR of the TS gene and exon 6 (2522) C/T of 
the SLC19A1 gene predict the survival of advanced NSCLC 
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patients treated with pemetrexed. However, a large scale 
clinical trial is required to validate these findings. 

Introduction 

Lung cancer is currently the most frequent cause of cancer 
mortality in the world and non-small-cell lung cancer (NSCLC) 
accounts for ~80% of all lung cancer cases (l).Moreover, the 
5-year survival of NSCLC is only 15% (2), as most patients are 
diagnosed at an advanced stage of the disease when surgery 
options are limited. Thus, chemotherapy has become the 
most important treatment option for treating these patients. 
Platinum-based chemotherapy is considered to be the founda- 
tion of treatment for patients with advanced NSCLC, with a 
median overall survival (OS) of 8-11 months, and therefore 
may prolong survival and improve the quality-of-life for these 
patients (3). However, current chemotherapeutic approaches 
have had limited success and therefore novel agents are 
urgently required to improve outcomes. Pemetrexed is a 
United States Food and Drug Administration (FDA) -approved 
therapy for advanced NSCLC (4) and a large phase III study has 
demonstrated that pemetrexed in combination with cisplatin 
provides similar efficacy and improved tolerability compared 
with gemcitabine in combination with cisplatin for advanced 
NSCLC patients (5). The study also indicated that pemetrexed 
provides similar improvements for patients with adenocarci- 
noma. Pemetrexed is widely used for first- and second-line 
and maintenance therapy in NSCLC (4-6). However, in China 
Pemetrexed is considered to be an expensive drug and is there- 
fore not widely used. Consequently, it is critical to identify 
patients who may benefit from pemetrexed therapy in order 
to reduce the side-effects and provide the most cost-effective 
approach. A number of clinical trials have attempted to iden- 
tify molecular biomarkers which may be associated with the 
response and toxicity of pemetrexed (7,8). The ultimate goal 
of this type of research is to create individualized treatment 
that provides an improved therapeutic profile by identifying 
the patients who are most likely to benefit from personalized 
therapy. 
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Pemetrexed is a multi-targeted, folate-based anti-metabolite 
that has been shown to have anticancer activity in other tumor 
types, including cervical (9), pancreatic (10), colorectal (11), 
breast (12) and gastric (13) cancers. Pemetrexed inhibits 
several enzymes in the thymidine and purine biosynthetic 
pathways, including thymidylate synthase (TS), dihydrofo- 
late reductase (DHFR), glycinamide ribonucleotide formyl 
transferase (GARFT) and aminoimidazole carboxamide 
ribonucleotide-formyl transferase (AICARFT). Pemetrexed is 
transported into the cell by the reduced folate carrier (RFC) 
system and is then rapidly and extensively polyglutamated by 
folylpolyglutamate synthase (FPGS). The drug is eliminated 
from the cell by y-glutamyl hydrolase (GGH). 

TS is a folate-dependent enzyme and the presence of a 
single nucleotide polymorphism (SNP) in the TS promoter has 
been shown to be an independent factor in metastatic colorectal 
cancer patients treated with 5-fluorouracil (5-FU) (14). 
Methylenetetrahydrofolate reductase (MTHFR) is a key enzyme 
in the folate metabolism pathway and MTHFR polymorphisms 
have been reported to correlate with the outcome of patients 
treated with methotrexate (15) and pemetrexed (16,17). The 
reduced folate carrier SLC19A1 is responsible for the transport 
of pemetrexed and therefore genetic variations in the gene may 
affect the transport of the drug (18). The associations between 
the gene polymorphisms and drug efficacy remain controver- 
sial. At present, few studies have investigated the correlation 
between pemetrexed pathway-associated gene polymophisms 
and drug efficacy. Therefore, the present study examined the 
association between the gene polymorphisms and therapeutic 
results of pemetrexed treatment in NSCLC patients. 

Patients and methods 

Patient selection. Chemotherapy-naive patients who were 
>18-years-old with histological or cytological evidence of 
measurable metastases or stage IIIB or IV NSCLC were 
enrolled in the present study. Other inclusion criteria included 
an Eastern Cooperative Oncology Group performance 
status (ECOG PS) <2 and normal organ function, including 
adequate hepatic, renal and hematological function. Patients 
who had brain metastases were eligible if they had a life 
expectancy > three months. Pregnant or lactating women 
were excluded from the study. All patients provided written 
informed consent prior to the initiation of treatment. Separate 
informed consent was obtained for the collection of blood 
samples for the study of single-nucleotide polymorphisms 
(SNPs). This study was approved by the ethics committee of 
the Second People's Hospital of Lianyungang (Jiangsu, China). 

Treatment . Pemetrexed (500 mg/m 2 ) was administered to 
patients as a 10 min intravenous (i.v.) infusion on day 1 of a 
21 day cycle with cisplatin (75 mg/m 2 , i.v.) administered on days 
1-3 every three weeks. All patients received daily oral supple- 
ments of folic acid (350-1000 pig) 5-7 days prior to the first day 
of administration and the supplementation was continued for 
three weeks after the treatment ended. Vitamin B 12 (1000 fig) 
was administered by intramuscular (i.m.) injection every three 
cycles of treatment. Dexamethasone (4 mg) was administered 
orally twice a day on the day before, day of and day after each 
dose of pemetrexed therapy to prevent allergic reactions. 



Treatment evaluation. All patients received complete clinical 
and laboratory evaluations prior to each cycle of treatment, 
including patient histories, physical examination, complete 
blood count, serum chemistry profile, electrocardiography 
(ECG) and radiological evaluation. Computed tomography 
(CT) and/or magnetic resonance imaging (MRI) were 
performed to evaluate the tumor response after every two cycles 
of therapy. The treatment effect was assessed according to the 
Response Evaluation Criteria in Solid Tumors (RECIST) (19) 
and categorized as complete response (CR), partial response 
(PR), stable disease (SD) or progressive disease (PD). 

DNA extraction. A 5 ml aliquot of blood was collected from 
each patient into an EDTA-coated tube before they underwent 
the first drug cycle. Genomic DNA was extracted from periph- 
eral blood cells using the Blood DNA kit (Omega Bio-Tek, 
Norcross, GA, USA) according to the manufacturer's instruc- 
tions. 

TS genotyping. DNA was amplified by polymerase chain 
reaction (PCR) using the following primers: 
5 -AACTGTGCTGCT GGCTTAGA-3' (forward) and 
5 -GTCTGTAAGGCGAG GAGGAC-3' (reverse). The TS 
gene promoter repeat and SNP (two or three repeats; G>C 
within the second repeat of the 3R allele) polymorphisms 
were directly detected by DNA sequencing using an ABI 
Prism 3730 DNA analyzer (Applied Biosystems, Carlsbad, 
CA, USA). The TS genotype was divided into low expression 
genotypes (2R/2R, 2R/3C and 3C/3C) and high expression 
genotypes (2R/3G, C/3G and 3G/3G) according to the SNP 
and TS genotype. 

MTHFR and SLC19A1 SNP. MTHFR and SLC19A1 SNPs 
were detected using a Taqman minor groove binder (MGB) 
probe-based assay. The PCRs for the MTHFR and SLC19A1 
genotypes were performed with a total volume of 5 jtl, 
including 2.5 jil Taqman Master Mix, 0.25 fA Taqman MGB 
primer probe, 1 pi\ DNA and 1.25 fi\ DEPC H 2 0. The PCR 
cycling parameters were as follows: 45 cycles of 95 °C for 
10 min, followed by 92°C for 15 sec and 60°C for 90 sec. 
Genotyping was performed using the ABI 7900 HT Sequence 
Detector (Applied Biosystems) and analyzed by the Allelic 
Discrimination program (Applied Biosystems). 

Statistical analysis. SPSS version 13.0 (SPSS, Inc., Chicago, 
IL, USA) was used for the data analyses. The primary end 
point was progression-free survival (PFS) and the secondary 
end points were objective response rate (RR) by RECIST 
and OS. PFS was measured from the time of randomiza- 
tion to disease progression or mortality from any cause, 
whichever occurred earlier. OS was measured from the time 
of randomization to the time of mortality from any cause. 
The Hardy-Weinberg equilibrium for gene polymorphisms 
was used to compare the observed distributions of genotype 
frequencies. Fisher's exact test was used to compare the treat- 
ment RRs between the genotypes. Kaplan-Meier curves and 
log-rank tests were also used to compare the OS and PFS 
distributions between the different SNP subgroups. P<0.05 
was considered to indicate a statistically significant differ- 
ence and all values were two-sided. 



Table I. Clinical characteristics. 



ONCOLOGY LETTERS 5: 1165-1170, 2013 

Table III. Response rate and gene polymorphisms. 



1167 



Characteristic 



Data 



Gene polymorphism 



CR + PR PD + SD P-value 



Gender, n (5) 

Male 

Female 
Age (years) 

Median age 

Range 

ECOG, n (%) 
0 
1 

2 

Stage of disease, n (%) 
IIIB 
IV 

Smoking status, n (%) 

Former/current smoker 

Never-smoker 
Number of cycles, n 

Total 

Median 

Range 



TS 



0.044 



23 (51.1) 


2R/3G, 3C/3G, 3G/3G 


7 


23 




22 (48.9) 


2R/2R, 2R/3C, 3C/3C 


8 


7 






MTHFR 






0.306 


63 


CC 


8 


9 




39-81 


CT 


5 


16 






TT 


2 


5 




12 (26.7) 


SLC19A1 






0.701 


24(53.3) 


CC 


3 


8 




9 (20.0) 


CT 


9 


14 






TT 


3 


8 





21 (46.7) 
24 (53.3) 

31 (68.9) 
14(31.1) 

199 
4 

1-6 



ECOG, Eastern Cooperative Oncology Group. 



CR, complete response; PR, partial response; PD, progressive dis- 
ease; TS, thymidylate synthase; MTHFR, methylenetetrahydrofolate 
reductase. 



Lianyungang hospital, which is affiliated with Bengbu Medical 
College. Two patients refused additional treatment after the 
first cycle administration due to grade 4 neutropenia. There 
were no significant differences in the baseline characteristics 
of the 45 patients and all the patients had adenocarcinoma. 
The median age was 63 years (range, 39-81 years). The char- 
acteristics of all the patients are shown in Table I. 



Table II. Distribution of TS, MTHFR and SLC19A1 genotypes. 



Genotype 


db SNP id 


N (%) 


TS 


rs45445694 




2R/2R 




0(0) 


2R/2C 




7 (15.6) 


3C/3C 




8 (17.8) 


2R/3G 




12 (26.7) 


3G/3G 




4 (8.9) 


3C/3G 




14 (31.1) 


MTHFR C677T 


rsl801133 




CC 




17 (37.8) 


CT 




21 (46.7) 


TT 




7 (15.6) 


SLC19A1 exon 6 (2522) C/T 


rsl051298 




CC 




11 (24.4) 


CT 




23 (51.1) 


TT 




11 (24.4) 



SNP, single nucleotide polymorphism; TS, thymidylate synthase; 
MTHFR, methylenetetrahydrofolate reductase. 



Results 

Characteristics. Between August 2009 and December 2011, a 
total of 47 patients were enrolled in the trial conducted at the 



Genotypes. The gene information was obtained from the 
HapMap database. Table II shows the genotype distribu- 
tions of TS, MTHFR and SLC19AL The distributions of the 
genetic variants followed the Hardy-Weinberg equilibrium 
(P>0.05). 

Association between gene polymorphisms and the treatment 
response to pemetrexed. There were 14 patients in the study 
who achieved a PR, eight with PD, 22 with SD and one who 
achieved a CR. The overall response rate (ORR; CR + PR) 
was 33.3% (15/45). All 45 patients had evaluable data for 
analysis. No significant associations were observed between 
the response rate and polymorphisms of MTHFR and 
SLC19A1. Patients with 2R/2R, 2R/3C or 3C/3C genotypes 
had a significantly higher RR than patients with 2R/3G, 
3C/3G or 3G/3G genotypes (53.3 vs. 23.3%, respectively; 
P=0.044; Table III). 

Survival analysis in patients with gene polymorphisms . 
According to the Kaplan-Meier survival analysis, patients 
with the TS 2R/2R, 2R/3C or 3C/3C genotypes had a signifi- 
cantly longer median PFS than patients with the 2R/3G, 
3C/3G or 3G/3G genotypes (6.8 vs. 3.8 months, respec- 
tively; log-rank x2 = 4.417; P=0.036; Fig. 1A). However, 
no significant difference in OS was observed between the 
two groups (P=0.638; Fig. 2A). The Kaplan-Meier survival 
analysis showed that PFS and OS did not differ for the three 
genotypes of MTHFR (PFS: CC, 5.6 months; CT, 3.8 months; 
TT, 5.8 months; log-rank x2 = 0.016; P=0.992; OS: 10.3 vs. 
10.6 vs. 8.1 months, respectively; P=0.739; Fig. IB and 2B). 
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Table IV. PFS and OS according to SNPs. 



Gene genotype 




PFS (months) 






OS (months) 




Median 


95% CI 


P-value 


Median 


95% CI 


P-value 


TS 






0.036 






0.638 


TD/on innp '3/~ , /'3/~ , 

zK/zK, ZK/3C, 3C/3C 


0.5 






1 n q 
1U.3 


o.j-lz. / 




iD/ifj irrar. ^nnr; 

ZK/JU, JWJO, ju/jU 


J .o 


] C\ 6, 6, 




1 0 1 
1U.1 


/ .o-lZ.o 




MTHFR 






0.992 






0.739 


CC 


5.6 


3.7-7.5 




10.3 


7.6-13.0 




CT 


3.8 


1.6-6.0 




10.6 


4.6-16.6 




TT 


5.8 


0.0-12.5 




8.1 


4.5-11.7 




SLC19A1 






0.137 






0.022 


CC 


4.9 


2.0-7.8 




12.2 


8.6-15.8 




CT 


5.8 


2.0-9.6 




8.9 


4.5-13.3 




TT 


5.5 


2.4-8.6 




7.3 


0.9-13.7 





PFS, progression-free survival; OS, overall survival; SNP, single nucleotide polymorphism; TS, thymidylate synthase; MTHFR, methylenetet- 
rahydrofolate reductase. 



A SnrvtriU fraction B Survival AmkHom C SurvtnlfttHiM 




ioo iooo hoo aooa •. m iooo iioo xao too noo no> aooo 

OS (months) OS (moatbs) OS (months) 

Figure 1. Kaplan- Meier curves estimating the association between the SNPs of the three genes and PFS of NSCLC patients treated with pemetrexed. 
(A) MTHFR, (B) SLC19A1 and (C) TS. SNPs, single nucleotide polymorphisms; PFS, progression-free survival; NSCLC, non-small cell lung cancer; TS, 
thymidylate synthase; MTHFR, methylenetetrahydrofolate reductase. 




Figure 2. Kaplan-Meier curves estimating the association between the SNPs of the three genes and the OS of NSCLC patients treated with pemetrexed. 
(A) MTHFR, (B) SLC19A1 and (C) TS. SNPs, single nucleotide polymorphisms; OS, overall survival; NSCLC, non-small cell lung cancer; TS, thymidylate 
synthase; MTHFR, methylenetetrahydrofolate reductase. 



No significant association was observed between the PFS of 
patients with the SLC19A1 CC genotype and patients with 
the CT and TT genotypes (log-rank x2 = 3.976; P=0.137; 
Fig. 1C). However, patients with the SLC19A1 CC genotype 



had a longer OS compared with the other two groups (12.2 
vs. 8.9 and 7.3 months, respectively; P=0.022; Fig. 2C). The 
PFS and OS summaries based on SNP analysis are shown in 
Table IV. 
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Discussion 

The present study was conducted to investigate the association 
of polymorphisms in pemetrexed-related genes with the clinical 
outcome in advanced NSCLC. The results suggest that the poly- 
morphisms in the TS 5'-untranslated region (UTR) may result in 
differences in the PFS of NSCLC patients treated with peme- 
trexed-based chemotherapy. The response rate was higher in 
patients with the TS 2R/2R, 2R/3C or 3C/3C genotypes. Patients 
with the SLC19A1 exon 6 (2522) CC genotype had improved 
OS when treated with pemetrexed compared with the patients 
without the genotype. However, a large-scale randomized 
prospective clinical trial is required to confirm these findings. 

The TS gene is located on chromosome 18pll.32 and the 
enzyme encoded by it is essential for cell proliferation as it 
catalyzes the methylation of deoxyuridine-5'-monophosphate 
(dUMP) to form deoxythymidine-5'-monophophate (dTMP). 
It is known that TS is a potential target for fluoropyrimidine 
and other drugs, including capecitabine, raltitrexed and peme- 
trexed. The promoter enhancer region of the TS gene is prone to 
polymorphisms with double or triple repeats of 28-bp tandem 
repeats in the 5'-UTR. Mandola et al reported that the 3R alleles 
have a G>C SNP at the 12th position of the second repeat (20). 
Based on this SNP and the TS genotype, the 2R/3G, 3C/3G 
and 3G/3G genotypes were considered to be high expression 
types and the 2R/2R, 2R/3C and 3C/3C genotype were consid- 
ered to be low expression types. However, no consensus was 
reached on the correlation between the TS polymorphisms and 
efficacy of 5-FU. Several studies have noted that a C SNP was 
associated with longer PFS in colorectal cancer patients treated 
with 5-FU (21,22). However, several other studies did not find 
any significant difference between TS SNPs and the clinical 
outcome of patients treated with 5-FU (23-25). 

Few studies have investigated the association between TS 
genotypes and the efficacy of pemetrexed. Smit et al showed 
that there was no correlation between the high and low TS 
expression genotypes and the clinical outcome of advanced 
NSCLC patients treated with pemetrexed (16,17). In the present 
study, the PFS was significantly longer in patients with the 
2R/2R, 2R/3C or 3C/3C genotypes compared with patients with 
the 2R/3G, 3C/3G or 3G/3G genotypes (P=0.036). However, no 
significant differences in OS were observed between the two 
groups (P=0.638). Further studies in larger populations are 
required to confirm these results, however, possible explana- 
tions include differences in the schedules of pemetrexed-based 
chemotherapy, types of lung cancer and ethnic populations. 

MTHFR is a central regulatory enzyme involved in folate 
metabolism and the cellular concentrations of 5,10-methylene- 
tetrahydrofolate are regulated by the activity of MTHFR (26). 
Moreover, optimal TS inhibition occurs when the cellular 
concentrations of 5,10-methylenetetrahydrofolate are 
elevated (27). The most common gene polymorphism linked 
with enzyme activity is MTHFR&ll C/T (28). At present, the 
association between MTHFR SNPs and the efficacy of 5-FU is 
inconclusive and few studies have investigated the association 
of SNPs in this gene with the efficacy of pemetrexed. Jakobsen 
et al reported that the RR was significantly higher in patients 
with the MTHFR TT genotype compared with two other 
genotype groups in metastatic colorectal cancer patients treated 
with 5-FU (29). However, several large-scale clinical trials 



did not identify a clear association (30-32). Moreover, a large 
meta-analysis indicated that MTHFR 677C/T is not a reliable 
predictor of the efficacy of FU-based therapy (33). Smit et al 
reported that patients with an MTHFR TT homozygous mutation 
had increased PFS compared with wild-type or heterozygous 
NSCLC patients treated with pemetrexed [7.9 months, (95% 
CI, 3.9-16 months) vs. 2.9 months, (95% CI, 2.8-3.4 months); 
P=0.03] in (16). However, Argiris et al reported contrasting 
results (17). In the present study, no correlations were observed 
between the MTHFR genotypes and efficacy of pemetrexed in 
advanced adenocarcinoma lung cancer. It is possible that the 
present results were due to the small sample size of the study or 
due to only one SNP in MTHFR being investigated. This may 
have prevented the identification of possible pharmacogenetic 
associations, since drugs may exert their anticancer effects 
through multistep, multigenic cascades (30). 

Pemetrexed is transported into cells by the RFC protein, 
which is encoded by the SLC19A1 gene. In previous studies, 
the SLC19A1 exon 6 (2522) CC genotype was observed to be 
associated with longer PFS and three-month progression-free 
status (progression-free at three months) in NSCLC patients 
treated with a combination of pemetrexed and bevacizumab 
chemotherapies (34). Adjei et al also demonstrated that a 
polymorphism in SLC19A1 may predict survival differences 
in pemetrexed-treated NSCLC patients (18). In the present 
study, patients with the SLC19A1 CC genotype had a longer 
median OS than patients with the other two genotypes (12.2 
vs. 8.9 and 7.3 months, respectively; log-rank x 2 = 2.957; 
P=0.022). The present results are similar to those of the two 
previously described trials. 

In summary, polymorphisms in the 5'-UTR of the TS 
gene and exon 6 (2522C/T) of the SLC19A1 gene appear to 
be potential prognostic factors for NSCLC patients treated 
with pemetrexed. Notably, germline polymorphisms have 
been reported to correlate with response and/or toxicity. The 
ultimate aim of this research is to aid the identification of the 
patients who are most likely to benefit from pemetrexed-based 
chemotherapy. Therefore, SNPs may be a useful tool for 
creating optimal individualized treatment regimens. 
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